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Abstract. We study the behavior of the solution of a generalized damped KdV equation Ut+Ux + 

Uxxx + u^Ux + = 0. We first state results on the local well-posedness. Then when p > 4, 

conditions on are given to prevent the blow-up of the solution. Finally, we numerically build 
such sequences of damping. 
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Introduction 

The Korteweg-de Vries (KdV) equation is a model of one-way propagation of small amplitude, 
long wave [KdV95] . It is written as 

Ui -\- Ux J- Uxxx T UUx — 0* 

In |BDKM96] . Bona et al. consider the initial- and periodic-boundary-value problem for the 
generalized Korteweg-de-Vries equation 

Ut + Uxxx + uPux = 0 

and study the effect of a dissipative term on the global well-posedness of the solution. Actually, 
they consider two different dissipative terms, a Burgers-type one —6uxx and a zeroth-order term 
au. For both these terms, they show that for p > 4, there exist critical values 5c and Uc such that 
if 5 > 5c or cr > fJc the solution is globally well-defined. However, the solution blows-up when the 
damping is too weak as for the KdV equation |MMn2| .The literature is full of work concerning the 
dampen KdV equation with p = 1 |ABS89[ ICR041 ICS13bl ICS13al IGhi88l IGhi941 IGouOOl IGR02j , 
but few are concerning more general nonlinearities. 

In our paper, we consider a more general damping term denoted by Our purpose is 

to find similar results as above, both theorically and numerically. So the KdV equation becomes 
a damped KdV (dKdV) equation and is written 

Ut Ux Uxxx T Ux ~\~ — 0. 

The damping operator ^^{u) works on the frequencies. It is defined by its Fourier symbol 

^^{u){i) := 7 ( 0 «( 0 - 
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Here u is the Fourier transform of u and 7 a strictly positive function chosen such that 

f u{x)^^{u)dfi{x) = [ > 0 . 

Jr jk 

We notice than the two cases studied in [BDKM^ are present with this damping by taking 
7(0 = 7 (^) = a respectively. 

The KdV equation has an infinite number of invariants such that the L^-norm. But, for the 
dKdV equation, the L^-norm decreases. Indeed, for all t € M, 


A 

dt 



u 


|2 

I7 


where the natural space of study is 


7/7 (M) 



s.t. 


[ 7(0 < +00j 

Jr J 


and the associated norm is 

\u\^ ■= ^_^7(0 

An other property of the KdV equation is that the solution can blow-up as soon as p > 4 . 
The blow-up is caracterized by lim l|rt||j 7 i = -|-oo. 

In this paper, we first establish the local well-posedness of the dKdV equation. Then we 
study the global well-posedness. More precisely, we focus on the behavior of the Fl^-norm with 
respect to p and we obtain conditions on 7 so there is no blow-up. Finally, we illustrate the 
results using some numerical simulations. We first find a constant damping ( 7 (.^) =constant) 
such that there is no blow-up and then the damping is weaken in such a way lim 7 (^) = 0 . 

l«|-5-+00 


1 Preliminary results 


Some results of injection concerning the space are given. 

Proposition 1.1. Assume < -t-00 then there exists a constant C > 0 such that ||u||oo < 

C\u\^, i.e., the injection H^{R) ^ is continuous. 

Proof. Let u S H.y(S.). We notice that 

u{x) = [ 

Jr 


Then 

< f\m\ = 

We assumed that 7 (^) > 0. Hence, the Cauchy-Schwarz inequality involves for all x € M : 


|tt(x)| < 


7(0 


7(01^(01' 


1 



□ 


Proposition 1.2. Let 7 and jd he such that for all ^ € M, 7 (^) > /3(^). We define 

m 


The continuous injection H.. 


p{N) := max . 

i>N 7(0 

FIs(M) is compact if and only if lim p{N) = 0. 

N ^+00 
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Proof. The condition is necessary. Indeed, if there exists a > 0 such that p{N) > a, ^N, then 
the norms |u|^ and |u|^ are equivalent, the injection cannot be compact. Let us prove now that 
the condition is sufficient. First, we have for u € H.y(R) : 

w\p= f m)m)\^< f 7(okop = K- 

JR JR ' 

This shows that the injection is continuous. Now we prove that the injection is compact. We 
use finite rank operators and we take the limit. Let In he the orthogonal operator on the 
polynomials of frequencies such that —N < I, < N. We have 


Inu = 






Thus 


\{id-iN)u\i= f 


< 


\i\>N 

f M 

'\i\>N 7(0 
|2 


7(0I^(0P 


<p(N)|u|; ^ 0 . 

^ A^^+oo 

Therefore Id is a compact operator and consequently the injection is compact. 


□ 


Proposition 1.3. Assume that u, v € H.y{ 

V^, rj we have 

< c - v) + . 


and there exists a constant C > 0 such that 


Then we have 


\uv\^ < C (\u\^ llullii + |u|..^ ||w||li) • 


Moreover if fj^ < +oo then II.y(R) is an algebra. 
Proof. Let u, v £ H.y(R). We have 


= f 7(0|S^(0P- 

JR 


We remind that rtu(0 = u* v{f). Using the inequality 

< C (^\J-f{^-ri) + , 

we obtain for all ry € M 

< C ( j \llii - r])\u{i - r])v{r])\dr] + f - v)v{r])\dri\ . 

\i/M / 

Hence 

\uv\‘^<c‘^ f (f - v)H^ - v)v{v)\dri + f - v)v{'n)\dT]] df,, 

t/M. \t/M. t/M / 

< 

< 


\ll{i - d)\HC - d)v{d)\dr]^ + \/lid)\u{C - rj)v{rj)\dr] 

\/7^l«l) * 1^ 




7(01^ 
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However, for f ^ and g ^ Lp', we have 

\\\f\*\9\\\h <Mh\\f\\h- 

Thus 

\uv\^ < C (iwl^ll'Ollii + I'yl^ll'wllii) • 

From DroDosition ll.il we know there exists a constant c > 0 such that HuHj;,! < c\u\^ if /j^ < 

+ 00 . Then, there exists (7 > 0 such that 

\uv\^ < 

□ 


2 Local well-posedness 


We study the following Cauchy problem : Vx G M, Vt > 0, 

f T T '^XXX T U^Ux T — 0, (1) 

j u(x,f = 0) = no(x). (2) 

The semi-group generated by the linear part is written as 

Stu-.= f 

JR 


In the rest of the section, f{u) denotes the non-linear part of the equation, i.e., f{u) = vPux- 
We first state a result of regularization. 


Lemma 2.1. Assume that s, r £ M^. Then there exists a constant Cr > 0, depending only on 
r, such that \/u £ H-yB[W) and Vt > 0 we have 


|5'tu| 


2 

-^s + r 



u 


2 

ryS • 


Proof. Let r £ M'*', 


u £ H.ys(M) and t > 0. Then we have 


I StU\^s+r 



2 

df 


< sup 
CgR 





But G M 


Thus 





a 

~p' 


l^iul 



u 


2 

7" ■ 


□ 


2 

Theorem 2.2. Assume that there exists r G]0, 2[ and for all ^ G M, 7 (^) > . We also assume 

that /g < -|-oo and there exists a constant C > 0 such that V^, r/ G M and s £ M"*" we have 

< c {yJiH - vY + ■ 

Then there exists a unique solution in ([—T, T], (E)) of the Cauchy problem (fT|)-(P|). 

Moreover, for all M > 0 with |uo|..^s < M and |uo|..j,s < M, there exists a constant Ci > 0 

such that the solution u and v, associated with the initial data uq and vq respectively, satisfy for 
2 

^ — \CoMp) 

W{-,t) - u(-,t)|.^, < Cl |uo - . 
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Proof. Thanks to Duhamel’s formula, ‘h(tt) is solution of the Cauchy problem, where 

$(u) = StUo - [ St-rf{u{T))dT. 

Jo 

Let show that u is the unique fixed-point of d>. We introduce the closed ball B(T) defined for 
T > 0 by 

B{T) := {[0,T];Hy{'R)) s.t. |n(t) - no(t)|^s < 3 |Mo|-y^} • 

We apply the Picard fixed-point theorem. We first show that (^B{T)^ C B{T). Let us take 
u S B{T) and show that <l>(u(t)) € B{T). We have 


< IS'tMol.yo + J \St-Tf{u{T))l 


On the one hand, we have 


\StUo\L = / 7(0^ StUo < / 7(0* k'ol < kolL • 


On the other hand, we apply Lemma l2.ll 


|5t_^/(u(r))|.^, = \St-rfiu{T))\^s-r+r 

\f{u{T))\^s-r 


< 


{t-ry. 


But 




{p + 1)2 7(^)' 


r7(0* 


UP+I 


df 


< 


w 


p+l 


(p +1 )^ 

2 

because 7(0 > O; and Hjs(R) beeing an algebra, we have 

\fiuiT))\ s-r < C'|n|^+^ . 


Consequently 


\^{u)\.ys < |uo|^, -h J 


dr 

0 (t-r )2 ^ 

* 1 


<koL*+C' sup [ - - :^dT 

te[o,T] ^ (t-r )2 

< |uo|.^* -k —sup • 

1-2 i6[0,T] ^ 2 


But u € B{T), then we have 

|u(t)|.^s - |uo|..y. < \u{t) - Uol^s < 3 |no|.^s • 


That involves 


and 


< 4|U0|.^. , 


sup < sup |u| « 

tG[0,T] ^ 2 \te[0,T] 


p+l 

< 4^+' |no|^+'. 
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We have <h(it(t)) € B{T) if the inequality 

\Hu{t)) - nol^. < 2 |uo|^. + ( 4 P +1 luol^t^) < 3 |no|, 

is true i.e. if 


0 < 2 < 


4P+^Cr\uoZs' 


Now let us show that <1> is a strictly contracting map. Let u, v a B(T), we prove that 

vte [o,r], 

sup \^{u{t)) - ^{v{t))\ s < k sup \u-v\ s 
ie[o,T] te[o,r] 

with k G [0,1[. As previously, we have 


Muit)) - <l>iv{t))\s = 


< 


/ St-r{f{u{T))-f{v{T)))dT 

lo 

r* cio 


/o {t — t) 2 


^P +1 _ yP+^ 


7 


Using the equality 

and the injection results, we obtain 


^p+i _ yP+i ^ 


= (u — v) ^ 

i+j=p 


UP+^ _ yP+^ 


<Ci\u — u| 


i+j=p 


< C 2 \U-V\^s '“'T- 

i+j=p 


U V 

u\ls lull. 


7“ 


< Cs |u - u| , |uo| 


p 

'■yii I '^U 


Then we have 


sup |4>(u(t)) - 4>(u(t))|.^, < C* |uo| 
te[o,r] 




'0 (t —r)2 


<C\uo\^sT^ 2 sup (|u-u| ,) 
te[o,r] ^ ^ 


The map is strictly contracting if 


1 T 

T^-2 < 


C\uoEs 


It remains to prove the continuity with respect to the initial data. Duhamel’s formula gives 
for t G [0,r], T^-i < 

\u - u|.^, < |uo - Uol.^, + 1/7) - f{v)\^s dr 

Wo - vo\^s + C'T^~^ ( W-v\y 


yi+j=p 


|uo - uol , + C"T^ 2 I Y I sup (|u-u|, 

\i+j=p j ie[ 0 T] ^ 


It involves 


|u - u| , < Cl |tto - UqL, . 


□ 
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Remark 2.3. Actually we can proove the local well-posedness for every 7 using a parabolic 
regularisation 

Ut Ux “ 1 “ Uxxx “ 1 “ A eUxx — 0 * 

Using lemma [KT\ with 7 (^) = the same computations as theorem \KE and taking the limit 
e —)■ 0 give the result jldrOtA I-B575]/ . 


3 Global well-posedness 

We work here under the hypothesis of the local theorem and study the global well-posedness of 
the damped KdV equation. We use here an energy method [BS751 IBS74| 

Theorem 3.1. If p < 4, for all 7 , the unique solution is global in time, valued in 77^(R). Else 
(p> A), there exists a constant 6 > 0 such that if > 0, G M then the unique solution is 
global in time, valued in 


Proof. Case p < 4: We begin by introducing N(u) and E{u), two invariants of the KdV equation 
without the damping term, which are the L^-norm and the energy. Their expressions are 


N{u) = / u^dx = ||m ||^2 , 

JR 


vF^'^dx = 7 \\ux\\\2 - 


(p+l)(p + 2 ) 


\P+2 


We first multiply (HD by tt and we integrate with respect to x. Then we have 

1 ^ 

2 dt 

Integrating with respect to time, we obtain 

rt 


■ / u^dx -|- / .^.y{u)udx = 0 . 
M. M 


u‘^dx + 2 [ ( [ AP.y{u)udx'\ dr = [ UqIx. 
: Jo Vjr / JR 


Which can also be written as 


V(tt)-I-2 J |u|^ dr = V(uo). 


We deduce from that expression that N{u) is a decreasing function and /q \u\idT is bounded 


independently of t by N{uo). Now, we multiply ([T]) by Uxx + and we integrate with respect 
to X. Then we have 


fL 

dt 


Ux 

(p + l)(p + 2 ) 


dx] — .^.y{ux)uxdx + / f- I dx = 0 . 

I Jr Jr \ P +1 / 


Integrating with respect to time, we obtain 


E{u) + f‘ W?, dr - l‘ dfj dr = E(u„). 


From this expression, we have 


/ Uxdx = E{u) -|- / 

< E{uo) + [ 

Jm 

< E{uo) + 


w 


,p+2 


{p+l){p + 2) 


U‘ 


P+2 




{p + l){p + 2) Jo 
1 „ „ 


Af.y{u) 


w 


,p+l 


P+1 


dx dr 


{p+l){p + 2) 


\\u \\^2 + sup ||n||^^ 

\0<T<t P+ 


- f ( f Af.y{u)udx] dr. 

+ 1 Jo \JR / 
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Using the inequality ||tt||^oo < 2 \\u\\j ^2 ||ii 3 ;|lL 2 and because /jj \^^{u)u\ = fjg^^y(u)u, 


u^dx < E{uq) + 


Since ||u ||^2 < ||iio|lL 2 


Then 


P 

22 


7 -TTTT-TTaN^IIt'Z ||“a:||r 2 T oup 

(p+l)(p + 2)" 0<r<t 


2+f II Ilf , / o- II 112 II 

"Ux\\r 2 + sup 22 \\u\\f 2 11^^: 


- - \ / “ I |2 

1^2 \\^X\\]^2 j j 


u^dx < Co + CiWuxWh + C 2 sup ||n3;||22- 

L 0<t<£ 


sup \\Ux\\i 2 -C sup \\Ux\\l 2 <Co. 

0<T<t 0<T<t 


(3) 


If there exists T > 0 such that lirn ||rtx|lL 2 = +00 then ||ua;|U 2 — C ||rta;||22 —t +00 since p < 4 

and this is impossible because of (l3|). Consequently, ||rta ;||^2 is bounded for all t and so is the 
//^-norm. 

Case p > 4: 

We estimate the L^-norm of Uxx- We multiply © with Uxxxx and we integrate with respect 
to X. Then we have 


r, 7, ( [ ^XX^"^ (d^)^xxxxdx - f vP UxVixxXxdX. 

2 at Vjr / Jm Jr 

Using two integrations by part, we have 

/ -^{u) '^XXXX dx — / (yUxx^'^xxdx. 

«/]R. «/M 

Let us work on the last term. Using integrations by part, we have 

5p 


(4) 


- [ vFuxUxxxxdx = —^ [ vF ^Uxulx-dx - p{p - 1) / “^uluxxdx. 


It follows that 


— / U‘^UxU 


P 

x^xxxxdx ^ TtII^IIoo ||oo II 7;^2 p{p ll'^^^lloo / \^x'^xx\dx. 

^ Jr 

But, from the Cauchy-Schwarz inequality 

/ ^ IIl^ IIl^ • 

Jr 


Then we have 


[ '^^^x'li’xxxxdx ^ Halloo ll^x II 00 II 7;^2 T p{p l)||^||oo ll^x ||oo ll^iT Hl^ ll^oia: Hl^ * 

Jr j 

Using the inequalty ||t£||^ < 2||r£||2,2||w 3 ;||^ 2 , we obtain 


/ UxUxxxxdx 


^{2\\u\\?2\\u 


■. P ~1 

3 1 \ 1 


1^2 II '^XX 11^2 


1^2 II ^31*^11^2 


1 X ^ 

i \ 2 


+p(p 1) ( 2||ri||^2112^2 ) II'^IIl^II '^xx IIl^ 


I l |2 
I'^ 3 : 3 ; Hi;, 2* 


'3^3? II 7;,2 • 





From ([3]), it leads to the inequality 


But 


Id,, „ 

- 7 ~ K^rr /■ 2 “h 

2dr 




< 0. 


(U^ 


— = 


^^{UxX)Uxx — ^Ua 


di 


= [ ( 7(0 - 

Jr 


The function O is increasing for its two arguments. We previously notice that t)||L 2 is an 
decreasing function with respect to the time. Then, if 7 (^) — O,\t=o > 0, \\uxxi-,'t)\\L^ does not 
increase for t > 0. Particularly, if 7 (^) > fl (||rio|| 2 , 2 , ||rtoa:a;||L 2 ) =• d, the semi-norm \\uxxi-,'t)\\L^ 
is bounded by its values at t = 0 . □ 


Remark 3.2. This result is also true on the torus T(0,L) where the operator is defined by 
its Fourier symbol 

^^{u){k) ■.= 'jkUk- 

Here Uk is the k—th Fourier eoeffieient of u and {'yk)kez o-re positive real numbers chosen sueh 
that 

/ u{x)J!f^{u)dfj.{x) = Vyfclufcp > 0. 


4 Numerical results 

In this part, we illustrate the theorem 13.11 numerically. Our purpose is first to find similar results 
as in |BDKM^ i.e. find a 7 ^ constant such that the solution does not blow up. Then build 
a sequence of 7 ^, still preventing the blow-up, such that lim 7 ^ = 0. Since dKdV is a low 

|A:|->-|-oo 

frequencies problem, we do not need to damp all the frequencies. 


4.1 Computation of the damping 

In order to find the suitable damping, one may use the dichotomy. We remind that our goal 
is to prevent the blow-up, i.e., avoid that lim ||rt||/fi = -|-oo. Let us begin finding a constant 

t^-l-CXD 

damping = yrt as weak as possible. We mean by weak that 7 has to be as lower as 

possible to prevents the blow up. Let 7 a respectively 7 e be the damping which prevents the 
explosion and which does not respectively. To initialize the dichotomy, we give a value to 7 and 
we determine the initial values of 7 a and 7 e. Then from these two initial values, we bring them 
closer by using dichotomy. The method is detailed in algorithm [1] and illustrated in Figures [1] 
and El 
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Algorithm 1 7 ^ and 7 e using dichotomy 


Require: 70 , e 

14 

end while 

1 

Initialisation of 7 : 7 = 70 

15 

7e = 7 

2 

Simulation with 7 fc = 7 

16 

7a = 27 

3 

if Explosion then 

17 

end if 

4 

while Explosion do 

18 

while 7 a - 7 e > e do 

5 

7 = 27 

19 

_ 7a+7e 

1 2 

6 

Simulation with 7 a: = 7 

20 

Simulation with 7 ^ = 

7 

end while 

21 

if Explosion then 

8 

7a = 7 

22 

7e = 7 

9 

7e = i 

23 

else 

10 

else 

24 

II 

11 

while Damping do 

25 

end if 

12 

’ 2 

26 

end while 

13 

Simulation with 7 a: = 7 




7fc 






/6 

If explosion with 70 


If damping with 70 



/e 

-> 


0 


7fc 



Figure 1: Initialization Figure 2: Dichotomy 

W extend the method to frequencies bands in order to build sequencies 7 ^ decreasing with 
respect to \k\ and tending to 0 when \k\ tends to the infinity. So we begin by defining the 
frequencies bands {Ni < N 2 < ■ ■ ■) and we proceed as previously but only on the frequencies 
1^1 > Ni. The method is described in algorithm [2] and illustrated in Figure [3] and 01 


Algorithm 2 7 ^ and 7 e on the band using dichotomy 


Require: 7 ^, N and Nb _ iter 

1 : Initialisation 7 = 7 a 
2: 7|A:|>Af = 0 
3: Simulation with 7 
4: if Damping then 
5: return 7a = 7 

6 : else 
7: 7 = 7a 

8 : while Damping do 

9: 7|fc|>Ar - — 2 — 

10 : Simulation with 7 

11 : end while 


12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 


7e = 7 

7a,|A:|>Ar = ‘^7\k\>N 

end if 

for i = 1 to Nh_iter do 

_ 7a,|fe|>JV+7e,|fe|>JV 

7|fc|>Af — 2 

Simulation with 7 
if Explosion then 

7e,|fc|>Ar = 7|fc|>Ar 

else 

7a,|A:|>Ar = 7|A:|>Ar 

end if 
end for 
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7fc 


I 


0 "i 


7 a obtained 
at step 2 

initial 7 a 
optimal 7 a 
optimal 7 e 

initial 7 e 


7fc 


0 


7 a obtained 
at step 3 
initial 7 a 

optimal 7 a 
optimal 7 e 
initial 7 e 


N2 


Figure 3: Initialization 


Figure 4: Find the damping 


4.2 Numerical scheme 

Numerous schemes were introduced in [CheSad]. Here we chose a Sanz-Serna scheme for the 
discretisation in time. In space, we use the FFT. Actually, the scheme is written, for all k, as 


1 + ~ + Ik) ) = 




ikAt 
p + l' 






(k). 


We find with a fixed-point method. In order to have a good look of the blow-up, we also 

use an adaptative time step. 


4.3 Simulations 

We consider the domain [—L, L] where L = 50. We take as initial datum a disturbed soliton, 
written as 


uo{x) = 1.01 X 


(P + l)(p + 2 )(c- 1 ) 


cosh 



p(c- 1 ) 


[x — ct — d) 


where p = 5, c = 1.5 and d = 0.2L. We discretise the space in 2^^ points. The Figure [S] shows 
the solution whithout damping, i.e., 7 ^ = 0, V/c. We observe that the L^-norm of Ux increases 
strongly and the solution tends to a wavefront (as in [BDKM96]). 
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10 

10 

10 


2.67906 


2.67906 


2.67906 


0 2 4 6 



t 


Figure 5: At left, solution at different times f = 0, 2, 4, 4.9925 and 5.3303. At right, Ff^-norm 
and L^-norm evolution without damping and a perturbed soliton as initial datum. Here p = 5. 


Using the methods introduced previously, we first find two optimal constant dampings 7 e = 
0.0025 and 7 a = 0.0027. As we can see in Figure El 7 e does not prevent the blow up. In the 
opposite in Figure [7] 7 a does. And we also notice that the two dampings are quite close. 




Figure 6 : At left, solution at different times f = 0 , 2 , 5, 10, 11 and 11.3253. At right, FF^-norm 
and L^-norm evolution with 7 ^ = 0.0025 and a perturbed soliton as initial datum. Here p = 5. 


12 








































solution at different times 



Figure 7: At left, solution at different times t = 0 , 2 , 5, 10, 15 and 20. At right, Ff^-norm and 
L^-norm evolution with 7 ^ = 0.0027 and a perturbed soliton as initial datum. Here p = 5. 

Considering more general sequences, particularly such that lim 7 ^. = 0. Using algorithm 

|fc|->-+oo 

El Figure [5] shows that the sequence ( 7 ^) as a frontier between the dampings which prevent 
the blow up and the other which do not. To illustrate this, we take two dampings written as 
gaussians. The first (denoted by 71 ) is build to be always above the sequence 7 a and the second 
(denoted by 72 ) to be always below. In Figures lUl and flUl we observe the damping 7 = 71 prevents 
the blow up. But if we take 7 = 72 , the solution blows-up. 



Figure 8 : Example of a build damping. Here the initial datum is the perturbed soliton. Here 
p = 5. 
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Figure 9: At left, solution at different times t = 0 , 2 , 5, 10, 15 and 20. At right, i/^-norm and 
L^-norm evolution with 7 = 7 i and a perturbed soliton as initial datum. Here p = 5. 



Figure 10: At left, solution at different times t = 0, 2, 5, 7 and 7.928. At right, F7^-norm and 
L^-norm evolution with 7 = 72 and a perturbed soliton as initial datum. Here p = 5. 

Conclusion 

We studied the behavior of the damped generalized KdV equation. If p < 4, the solution does not 
blow-up whereas if p > 4, it can. To prevent the blow-up, the term 7 defining the damping has 
to be large enough. In particular, we build a sequence of 7 which vanishes for high frequencies. 
This frequential approach for the damping seems useful for low frequencies problem. 

Acknoledgments. I would like to thank my thesis supervisors, Jean-Paul Chehab and Youcef 
Mammeri, for their help and comments. 
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